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We argue that the Lagrangian for gravity should remain bounded at large curvature, and interpo- 
late between the weak-field tested Einstein-Hilbert Lagrangian £eh = R/16nG and a pure cosmo- 
logical constant for large R with the curvature-saturated ansatz C cs — £eh/v1 + l 4 R 2 , where I is a 
length parameter expected to be a few orders of magnitude above the Planck length. The curvature- 
dependent effective gravitational constant defined by dC/dR — 1/1Q-kG c s is G c g = Gy/l + l A R 2 , 
and tends to infinity for large R, in contrast to most other approaches where G c s — > 0. The theory 
possesses neither ghosts nor tachyons, but it fails to be linearization stable. In a curvature saturated 
cosmology, the coordinates with ds 2 = a 2 [da 2 / B(a) — dx 2 — dy 2 — dz 2 ^ are most convenient since 
the curvature scalar becomes a linear function of B(a). Cosmological solutions with a singularity 
of type R — > ±oo are possible which have a bounded energy-momentum tensor everywhere; such a 
behaviour is excluded in Einstein's theory. In synchronized time, the metric is given by 

ds 2 =dt 2 -t 6/5 (dx 2 + dy 2 + dz 2 ). (0.1) 

On the technical side we show that two different conformal transformations make £ cs asymptotically 
equivalent to the Gurovich-ansatz C = |i?| 4,/3 on the one hand, and to Einstein's theory with a 
minimally coupled scalar field with self-interaction on the other. 
PACS 98.80; 04.50 



I. INTRODUCTION 

According to an old idea by Sakharov [jlj, the gravitational properties of spacetime are caused by the bending 
stiffness of all quantum fields in a spacetime of scalar curvature R. This idea of induced gravity has inspired many 
subsequent theories of gravitation, from Adler's [Q proposal to consider Einstein gravity as a symmetry breaking effect 
in quantum field theory to the modern induced gravity derived from string fluctuations Whatever the precise 
mechanism, any induced gravity will lead to a Lagrangian which is bounded at large R, and may also go to zero. The 
latter case would be analogous to the elastic stiffness of solids, which is constant for small distortions, but vanishes 
after the solid cracks. 

In this paper we investigate the physical consequences of a simple Lagrangian which goes to a constant at large R, 
thus interpolating between the Einstein-Hilbert Lagrangian for small R and a pure cosmological constant for large R. 
This Lagrangian will be referred to as curvature-saturated and reads 

1 /? 

£cs=— , (1-1) 

The length parameter I may range from an order of the Planck length lp or a few orders of magnitude larger than lp. 
Applying standard methods and those of Refs. @-@], we shall derive the cosmological consequences of the saturation 
and compare our ansatz with others. 
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One of the motivations for a renewed interest in a more detailed consideration of cosmology with non-linear curvature 
terms comes from M-theory, see Ref. |1 "Brane new world" . In 0] a conformal anomaly is considered, which turns 
out to have analogous consequences as Starobinsky's anomaly-driven inflation with R- and i? 2 -terms, see e.g. Refs. 
pOf for the older results. Ref. Q] contains the latest results concerning the effective A-term in such models. 

Our own direct motivation to tackle the model discussed below was as follows: We tried to make the analogy 
proposed in 0] more closer than done by others; the analogy with solid state physics is this one: For small forces, the 
resistance to bending is proportional to this force, but after a certain threshold - defined by cracking the solid - the 
resistance vanishes. 



A similar line of reasoning was deduced in Ref. |12|: There the finite-size effects from the closed Friedmann universe 
to the quantum states of fields have been calculated. Instead of continuous distribution of the energy levels of the 
quantum fields, one has a discrete spectrum. Qualitatively, the result is: If the radius a of the spatial part of 
spacetime shrinks close to zero, which is almost the same as very large R, then the spacings between the energy levels 
become larger and larger, and after a certain threshold, all fields will be in the ground state. This behaviour shall 
be represented by an effective action. The concrete form of the corresponding effective Lagrangian is not yet fully 
determined (that shall be the topic of later work) , but preliminarily we found out that the behaviour for large R will 
quite probably be of a Lagrangian bounded by a special effective A; so we have chosen one of the easiest analytic 
functions possessing this large- R behaviour together with the correct weak-field shape. 

The paper is organized as follows: In Sec. II we calculate the consequences of the effective Lagrangian C cs . 

In Sec. Ill we investigate the consequences of the i?-dependence of the effective gravitational constant defined by 

I dC 



167rG off dW 

which is 



(1.2) 



G cff = GVl + l 4 R 2 (1.3) 

for C — C cs and tends to infinity as R — > ±oo. 

Then we apply two different conformal transformations to £ cs . One of them, presented in Sec. IV, makes C cs 
asymptotically equivalent to the Gurovich-ansatz jl3), Q 

The other transformation, by the Bicknell theorem given in Sec. V, establishes a conformal relation to Einstein's 
theory, with a minimally coupled scalar field. In the literature, see fl5| ] and the references cited there, only the second 
of these conformal transformations has so far been used. The physical consequences of these three theories are, of 
course, quite different since the metrics are not related to each other by coordinate transformations. 

Our approach differs fundamentally from that derived from the limiting curvature hypothesis (LCH) in Refs. [ jl6| , 
where the gravitational Lagrangian reads 

£ = ^ + ^(Vi-#7A 2 -i) (1.5) 

whose derivative with respect to R diverges for R — > A. This divergence was supposed to prevent a curvature 
singularity, a purpose not completely reached by the model presented in the first of Refs. [ |l6| because other curvature 
invariants may still diverge. (Let us note for completeness: In the second of Refs. Jig ], a more detailed version of 
the LCH is presented which covers also the bounding of the other curvature invariants; it is restricted to isotropic 
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cosmological models. For more general space-times one faces the problem that sometimes a curvature singularity 
exists, but all polynomial curvature invariants remain bounded there.) 



In contrast to Eq. (1.5), our model favors high curvature values. 



It turns out that the use of synchronized or conformal time is not optimal for our problem. We therefore use a new 
time coordinate which we call curvature time for the spatially flat Friedmann model. The general properties of this 
coordinate choice are described in Sec. VI. 

In Sec. VII we study the consequences of curvature-saturation for some cosmological models using the coordinates 
of Sec. VI. In Sec. VIII, finally, we summarize our results and compare with the related papers jL7j to p8| . 

II. FIELD EQUATIONS OF CURVATURE-SATURATED GRAVITY 



The curvature-saturated Lagrangian (1.1) interpolates between the Einstein-Hilbert Lagrangian 

R 



£eh = 



16ttG" 



which is experimentally confirmed at weak fields, and a pure cosmological constant at strong fields 

C CS (R) 



1/16ttGZ 2 




*► R 



-1/lQ-wGl 2 



(2.1) 



FIG. 1. Curvature-Saturated Lagrangian as a function of the curvature scalar 7?. 

±1 



16ttGZ 2 ' 

The R dependence is plotted in Fig. |l|. 

The usual gravitational constant is obtained from the derivative of the Einstein-Hilbert Lagrangian: 



(2.2) 



dC 



EH 



(2 1) 

16ttG dR ' K ' ' 

From our curvature-saturated Lagrangian ([hi]) we obtain, with this derivative, the effective gravitational constant 



(1.3). The definition (2.3) is motivated as follows: If one considers the Newtonian limit for a general Lagrangian 
C(R) which may contain a nonvanishing cosmological constant, the potential between two point masses contains a 
Newtonian 1/r-part plus a Yukawa-like part exp(— r/ry) stemming from the nonlinearities of the Lagrangian; the 
details are given in the Appendix. At distances much larger than ry , but much smaller than 1 / vR, only the 1 / r-term 
survives, and the coupling strength of the 1/r-term is given by the effective gravitational constant G c ff. For a recent 
version to deduce such weak- field expressions, see Ref. [p4|. 



For a general Lagrangian L(R) such as (LI), the calculation of the field equation is somewhat tedious, since the 
Palatini formalism which simplifies the calculation in Einstein's theory is no longer applicable. Recall that in this, 
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metric and the affine connection are varied independently, the latter being identified with the Christoffel symbol only 
at the end. 

Here the following indirect procedure leads rather efficiently to the correct field equations. Let 



c' = — c" d2c 



(2.4) 



dR' ~ dR 2 ' 

and form the covariant energy-momentum tensor of the gravitational field which is given by the variational derivative 
of C with respect to the metric g ab : 



2 SCy^g 



V~9 og ab 

where g denotes the determinant of g ab . For dimensional reasons, & ab has the following structure 

©ab = aC'R ab + /3C'Rg ab + ^Cg ab + SO£'g ab + e£' b 



(2.5) 



(2.6) 



with the 5 real constants a . . . e. These constants can be uniquely determined up to one overall constant factor by the 
covariant conservation law 



Q ab . h = 0. 



(2.7) 



The overall factor is fixed by the Einstein limit I — > of the theory, where Q ab = {Rab — \Rg a b) /8nG. In this way we 
derive the following form of the covariantly conserved energy-momentum tensor of the gravitational field 

1 



&ab = 



16ttG 



{2£'R ab - Cg ab + 2UC'g ab - 2C'. ab 



(2.8) 



The calculation is straightforward, if one is careful to distinguish between (□£'). a and □(£! a ), which differ by a 
multiple of the curvature scalar. 

Inserting our curvature-saturated Lagrangian ( |1.1| ) into (2.4) and omitting the subscript, we have 

dC 1 



c = | (i + ; 4 i? 2 )- 1/2 



and find from ( 



Rab 



Rg a b 



8irG I (1 + 24^2)3/2 2 (1 + Z 4 i? 2 ) 



1/2 



{l + l^R 2 ) 



3/2 



(l + l 4 R 2 



,3/2 



(2.9) 



(2.10) 



- ;ab 



Setting I — reduces this to l/16nG times the Einstein tensor. The trace of (2.1C) is 

R + 2l 4 R 3 



<~> a — 



1 

8^G 



{1 + l^R 2 ) 



3/2 



3D 



1 



(1 + Z 4 i? 2 ) 3/2 



(2.11) 



According to Einstein's equation, Q ab has to be equal to the energy momentum tensor of the matter T ab , i.e., T ab = Q ab 



Equation (2.11) implies that in the vacuum, the only constant curvature scalar is R — 0, such that this model does 
not possess a de Sitter solution. Further, we can see from Eq. ( p. 10 ), that a curvature singularity does not necessarily 
imply a divergence of energy-momentum, but may be compensated by the infinity of G c g. 



III. EFFECTIVE GRAVITATIONAL CONSTANT AND WEAK-FIELD BEHAVIOR 



Let us compare the effective gravitational constant G c g of our curvature-saturated model with those of other models 



discussed in the literature. From (1.3) we see that G c g has the weak- field expansion 
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G e ff = G 1 + -Z 4 i? 



and the strong-field expansion 



The full i?-behavior is plotted in Fig. 



Geff = Gl 6 \R\ 3 I 1 



2l 4 R 2 



(3.1) 



(3.2) 




FIG. 2. Effective gravitational constant as a function of the curvature scalar. 



The weak-field expansion of C cs is given by 



r — 



R 



R 



OO 



16ttGV1 + l 4 R 2 16ttG ^ 



b k R 



2k+l 



(3.3) 



with real coefficients bk, where bi = —l 4 /32irG. 

As one can see, the quadratic term is absent, so that the linearized field equation coincides with the linearized 
Einstein equation. Thus we encounter neither ghosts nor tachyons; for details see Appendix B. 

There is, however, a price to pay for it. The theory has lost linearization stability of the solutions. This latter 
property has the following consequences: If one performs a weak-field expansion 



9ij = >},, + 



e 9i 



(3.4) 



around flat spacetime to solve the field equation, one has to use the terms up to the order m — 2 to get the complete 
weak-field part of the set of solutions. With this peculiarity, we obtain a well-posed Cauchy problem for the gravity 
theory following from the Lagrangian C cs . 

Let us now compare our theory with others available in the literature. Let 



£*ct,n (-^) 



R 



16ttG 



aR n 



(3.5) 



with some number n > 1 and constant a / 0. In analogy with Eq. (1.2) we calculate the effective gravitational 
constant from 



1 



167rG c ff 



dR 



1 



16ttG 



anR 11 " 1 



(3.6) 
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such that 



G 



(3.7) 



1 + 167^cmi? n ~ 1 G , 

i.e., G e g — > as R — > ±oo. For n — 2, more exactly: for all even natural numbers n, we meet an additional peculiarity 
that G e s can diverge for finite values of R already. Such values of R = i? cr it are called critical B. For n = 2 we get 



Rc 



1 



(3.8) 



32a7rG 

and this is the region where G e ff changes its sign, as shown in Figures 3 and 4. At critical values of the curvature 
scalar, the Cauchy problem fails to be a well-posed one. 




FIG. 3. Effective gravitational constant G e g for C a ,a with a > as a function of R. 




G 



IV. CONFORMAL DUALITY 



In Rcf. [gj, a duality transformation relating between different types of nonlinear Lagrangians has been found. In 
the present notation it implies the following relation. Let 



rii 

9ab — >~ 9ab 



(4.1) 



be the conformally transformed metric with £' ^ 0, which is fulfilled by our Lagrangian (1.1). Then the conformally 
transformed curvature scalar equals 



R = 



3R 4C 



and the associated Lagrangian is 



£ 



2R 3£ 
£73 _ -jjjl- 



(4.2) 



(4.3) 



We easily verify that £' £' = 1. Then one can prove that g a b solves the vacuum field equation following from £(R) if 
and only if g a f, of Eq. (IO) solves the corresponding equation for £(R) of Eq. (O). 



Example: For £ = R 



k+l 



find, up to an inessential constant factor, £ = R k+1 with k = 1/(2 — 1/k), such that 



for a purely quadratic theory with £ = R 2 , also £ = R 2 . For our curvature-saturated model £ — > const, we should 
expect a behavior with k — ► — 1, i.e., k — > 1/3, this leads to £ ~ i? 4 / 3 , which is the Gurovich-model jL3|, cf. Eq. (|l.-j|). 

Let us study this in more detail. To simplify the expressions we use, in this subsection only, reduced units with 
\QttG = 1 to best exhibit the fixed point / = of this transformation making it an identity transformation if applied 



to Einstein's theory where k = 1. In the present units, Eqs. (2.9) have to be multiplied by 2 and become 

£ = R(1 + l 4 R 2 y 1/2 , £' = ^§ = (1 + l 4 R 2 ) ~ 3/2 . (4.4) 



Inserting these into (|4.l|)- ([4.3[) , we obtain 



9ab 



9ab 



(l + l 4 R 2 ) 3 



and 



For small R we have 



and for large \R\ 



r = -R{i + rR^) A {i~±rR z ). 



R 



-R(l- PR 2 + ...), 



R = 4l 16 R 9 | 1 
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4Z 4 i? 2 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



The inverse function R(R) of (4.6) is not expressible in closed form, but its small- and large-curvature expansion can 
be calculated from (4/7) and (L; 



R= -R{l + TR Z + ...) 



R 



R 
ills 



1/9 



11 

36F 



4/16 \ 2 /9 

R 



(4.9) 



^From Eq. (4.3) we see that 
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£= -i?(l + Z 4 i? 2 ) 9 / 2 (l-3Z 4 i? 2 ) 
where R(R) has to be inserted. For large R we use the right-hand equation in (E 



£ = 3/ 



22 



R 

4/16 



4/3 



51 f Al 



6Z 4 V R 



16 \ 2 /9 



(4.10) 

and obtain the limiting behavior 

(4.11) 



V. BICKNELL'S THEOREM 



Bicknell's theore J25|] , in the form described in Ref. |4[], relates Lagrangians of the type (2.9) to Einstein's theory 
coupled minimally to a scalar field <f> with a certain interaction potential V((j)). This Lagrangian is given by 



£eh + \<t>>iV - v{4>) 



(5.1) 



The relation of V((j>) with £(i?) is expressed most simply by defining a field with a different normalization ip = -\/2/3 (f>, 
in terms of which the potential V(<fi) — V{ip) reads 



V(i/j) = £(i?)e 



-2V- _ ^ e -V- 



with i? being the inverse function of 



V> = ln[2£'(i?)]. 



The metric in the transformed Lagrangian (5.1) is 



9ab = e^gab- 



For our particular Lagrangian (2.E) we have from (5.3) 



^ = --ln(l + Z 4 i? 2 ). 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



Now we restrict our attention to the range R > where ijj < 0; the other sign can be treated analogously. Then (5.5) 
is inverted to 



such that (5.2) becomes 



R= iry/e-W/3-i, 
V 



V(iP) = -^(e^ 5 ^/ 3 - e-^y/e-W/* - 1. 



(5.6) 



(5.7) 



In the range under consideration, this is a positive and monotonously increasing function of —ip (see Fig. ^), with the 
large-$ behavior 

1 



-2ip 



2/ 2 



(5.8) 



This is the typical exponential potential for power-law inflation. As mentioned at the end of Section II, no exact de 
Sitter inflation exists. For ip -> 0, also V(tp) -> like 4 V / 273V' 3/2 . 
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FIG. 5. Potential V(ip) associated with curvature-saturated action via Bicknell's theorem. 



If V(ip) has a quadratic minimum at some i/jq with positive value Vo = V(ipo), then there exists a stable de Sitter 
inflationary phase. As a pleasant feature, the potential V{ip) has no maximum which have given rise to tachyons. 



From Eq. (5.5) one can see that for weak fields, tp ~ R 2 , whereas a R + i? 2 -theory has rp ~ R. In other words: 
In our model it is a better approximation to assume the conformal factor to be approximately constant for weak 



fields then in R + i? 2 -theories, since at the level keeping only terms linear in R the two metrics g a b and g a b in (5.4) 
coincide. 



VI. FRIEDMANN MODELS IN CURVATURE TIME 



The expanding spatially flat Friedmann model may be parametrized with the help of curvature time a > as 
follows: 

da 2 



ds z = a z 



B(a) 



dx — dy — dz z 



(6.1) 



(6.2) 



where B(a) is an arbitrary positive function determining R as 

R = --—, 

a 3 da ' 

depending only on the first derivative of B{a). This is a special feature of ( |6.l| ) since, in general, the curvature scalar 
depends on the second derivative of the metric components. Note also the linear dependence of R on B' = dB/da, 
in contrast to the usual nonlinear dependence of the curvature scalar on the first derivative of the metric coefficients. 

Let us recall some facts on Friedmann models in curvature time and exhibit the corresponding transformation to 
synchronized time. 



A. From curvature time to synchronized time 



The spatially flat Friedmann model in synchronized time has the metric 

ds 2 = dt 2 - a 2 (t)(dx 2 + dy 2 +dz 2 ). 



Metric (3.1) goes over to metric (6.3) via 



dt 



a da 

7W) 



(6.3) 



(6.4) 



such that 







The inverse function a(t) provides us with the desired transformation. 



B. From synchronized time to curvature time 

Consider a(t) in an expanding model with 

d^>0. (6.6) 

Then we can invert a(t) to t(a), and have 

B(a) = a 2 [a{t{a)} 2 . (6.7) 

^From this relation we understand why R depends on the first derivative of B only: B itself contains a derivative of 
a, and R is known to contain up to second order derivatives of ait). 

C. Examples 

Let a(t) = t n , i.e. t = a 1 /", a(t) = nt n -\ a(t(a)) = na 1 " 1 /". Then Eq. (^) yields 

B(a) =n 2 a A - 2/n . (6.8) 
Let further a(t) = e Ht , H = const. > 0, a = Ha. Then 

B(a) = H 2 a 4 . (6.9) 



Obviously, Eq. ( p.9j ) is a limiting form of Eq. (|6.8| ) for n — * oo. Equation ( J6 . If ) with B(a) from ( |6 .9|) represents a 
vacuum solution of Einstein's theory with A-term where A = 3-ff 2 , namely the de Sitter spacetimc. 
Let us also give some examples for the direct use the curvature time: 



1.) From Eq. ( |6.2D we see that R = implies B = const., corresponding to n = | in Eq. (|6.§|), i.e., a = t 1 ! 2 
synchronized time. This is the usual Friedmann radiation model. 



2.) Also from Eq. fl6.2| ), a constant R ^ implies £? = Ci + C2& 4 with constants Cj and C2, C2 7^ 0. 
For Ci = 0, C2 = H 2 , this represents the de Sitter spacetime Eq. 



3.) The dust-model in synchronized coordinates is given by a — t 2 / 3 , i.e., with Eq. (|6 .8|) we get 



B(a) = ta, (6.10) 



such that B' = const. Together with Eq. (6.2), this leads to 

i?a 3 = const, (6.11) 

ensuring mass conservation, because R is proportional to the mass density, and the pressure is negligible for 
dust. 
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D. The variational derivative 



For the metric (6.1) we have 



V-9 = y/-dct gij = -j=. 
The Lagrangian for Einstein's theory with A-term reads 

C = (i? + 2A)V^9 ■ 

With (fj) and we get from j6l^ ) 



The vanishing of the variational derivative 



6C dC f dC\ 



SB dB \ dB 







(6.12) 



(6.13) 



(6.14) 



(6.15) 



gives B = H 2 a 4 with A = 3if 2 , i.e., the usual de Sitter spacetime. No integration is necessary, since the derivative of 
B cancels. Intermediate expressions are 



dC 
dB 



2A- 



1 



£-3/2 



(6.16) 



dC 
~dB' 



-iaB- 1 ' 2 



(6.17) 



E. Remaining coordinate-freedom 



Translations in t do not change the form of the metric (|6.3|) . This freedom is related to the fact that the integration 



constant in the integral (3.5) remains undetermined; this coordinate freedom has no analog in the metric in curvature 
time Eq. (pjj). 



The metric (5.1) has the following property: It remains unchanged under multiplication of a and B by the same 



positive constant. Such a constant factor appears if we multiply the spatial coordinates by a constant factor. In 
synchronized coordinates this property means that not a itself, but only the Hubble parameter 



H(t) 



(6.18) 



has an invariant meaning. By the same token, not B(a) itself, but only B{a)/a A has an invariant meaning. In fact, 



from Eq. (6.7) we see that 



B 



H 



(6.19) 



VII. COSMOLOGICAL SOLUTIONS 



Here we recall some formulas of Ref. and present some new results for the curvature-saturated Lagrangian. 
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A. Solutions for Lagrangian R m 

For the Lagrangian C — R m , we obtain the following exact solutions for a closed Friedmann universe: 



t 2 

ds 2 = dt 2 - —: — 7 da?,s, (7.1) 

2m 2 -2m-l (+) ' 



where da 2 + ^ is the metric of the unit 3-sphere. 



Analogously, for the open model 



t 2 

ds 2 =dt 2 -- —5 -daf +) . (7.2) 

2m - 2m 2 + 1 (+) v ; 

Of course, both expressions are valid for positive denominators only. 

For the spatially flat Friedmann model, it proves useful to employ the cosmic scale factor a itself as a time-like 
coordinate. 

ds 2 = a 2 [Q 2 (a)da 2 - dx 2 - dy 2 ~ dz 2 ] . (7.3) 

This coordinate is meaningful as long as the Hubble parameter is different from zero, so that we cover only time 
intervals where the universe is either expanding or contracting. Possibly existing maxima or minima of the cosmic 
scale factor as seen in synchronized time can, however, been dealt by a suitable limiting process and patching. The 
curvature scalar reads now 

a 3 Q 3 da ' ^ 
and to reduce the order of the field equation it proves useful to define 

P(a) = (7.5) 
da 

Then the field equation is fulfilled if 

dP P 

= m(m - 1)— + (m - 1)(1 - 2m)P 2 + m(4 - 3m) — . (7.6) 
da a 

Therefore, the spatially flat Friedmann models can be solved in closed form, but not always in synchronized coordi- 
nates. 



B. Solutions for Lagrangian £ cs 

In the context of our curvature-saturated model, we shall restrict ourselves to the expanding spatially flat Friedmann 
model. The field equation written in synchronized or conformal time — the two most often used time coordinates used 
for this purpose — have the disadvantage that the number of terms is quite large, and that even in the simplest case 
C = \R 2 we cannot give closed- form solutions, apart from the trivial solutions R = having the same geometry as 
the radiation universe (a = y/i in synchronized time t) and the de Sitter universe (a = e* in synchronized time t). So, 
we prefer to work in the less popular coordinates ( |7.3| ). In principle, the field equation should be of fourth order, but 
we shall reduce it to second order. 

To find the field equation for a spatially flat Friedmann model with our Lagrangian, it is useful to consider first 



a general nonlinear Lagrangian and specialize to £ cs afterwards. To simplify (7.4), we define instead of Q(a) the 
function B(a) = Q(a)~ 2 > as a new dependent function. Then (|7.3|) reads 
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ds 2 



" da 2 

Ma) 



— dx — dy 



dz 2 



and (7.4) goes over to 



R=-- — 

a 3 da 



(7.7) 



(7.; 



Thus, B itself does not appear explicitly, and only first, and not second derivatives are present. The geometric origin 
of this property is the same as in Schwarzschild coordinates — one integration constant is lost in the definition of the 
coordinates, and this makes curvature depend only on the first derivative of the metric. 



From the 10 vacuum field equations (2.10) only the 00-component is essential; it is the constraint equation, therefore 



it has one order less than the full field equation, but if the constraint is fulfilled always, then all other components 
are fulfilled, too. Together with Eq. (^J) we should now expect that the fourth order field equation ( |2.10|) can be 
reduced to one single second order equation for B(a), where hopefully, B itself no more appears. 



The equation 0qo = is via (2.4) and ( |2.8| ) equivalent to 



dB 



= 3£' ( 2B- a— - a 4 £-l8aB£"— — — 
da J da \ a A da 



1 dB 



(7.9) 



which is much simpler than the analogous equation in synchronous time, as observed here for the first time. 

Before we insert our Lagrangian £ cs into (|7.9|) , let us cross check its validity by solving known problems: If C" 
vanishes identically, then C! is a constant, and we return to Einstein's theory. The case B = const, gives the radiation 
universe, while B = a 4 is the exact de Sitter solution. For the Lagrangian £ = ^R 2 with £ = R and C" — 1, and 



Eq. (7.9) reduces to 



= aB 2 - AaBB + 8BB, 



(7.10) 



where a dot denotes differentiation with respect to a. Again, B = a 4 is the exact de Sitter solution. Defining (3 = lnS 
and z = a$, Eq. ( 7.10| ) goes over in 



4ai = 3z(4 - z). 



With a = In a we arrive at 



^ dz 
da 



3z(4 



(7.11) 



(7.12) 



which can be solved in closed form. Qualitatively it is clear that z = 4, i.e., the de Sitter solution, represents an 



attractor. Solving Eq. (7.12) we obtain in the region < z < 4: 



z = 2 



2 tanh ( —a 
2 



showing explicitly that z — > 4 for a — ► oo. The metric can be calculated from 

i 3 -r 







i 



,3 



1 



(7.13) 



(7.14) 



using the identity 



1 This behavior is known already from the Friedmann equation in General Relativity: Energy density is proportional to the 
square of the Hubble parameter which contains only a first derivative. 
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tanh In x 



1 



1 



(7.15) 



After these preparations we are ready to deal with our Lagrangian £ cs . We insert C and C from Eq. (4.4), and 

C" = -3Z 4 i?(l + l A R 2 )- b ' 2 (7.16) 



into Eq. (7.9) and obtain, after setting I = 1, the simple expression 



da 



9B 2 )(2B - aB) + B(a 6 + 9B 2 ) 2 . (7.17) 
In these coordinates, the flat Minkowski spacetime does not exist, and the radiation universe R = is not a 



solution. This is why B = const, yields no solution to Eq. (7.17). Also, as was known from the beginning: the de 
Sitter spacetime B = a is not an exact solution here. However, in the nearby-region where the Lagrangian is well 
approximated by a quadratic function in R with a nonvanishing linear term, the behavior of the solutions is quite 
similar to that of R + i? 2 -models, where no exact de Sitter solution exists, but a quasi de Sitter solution represents a 
transient attractor with sufficient long duration to solve the known cosmological problems. These calculations have 
been presented at different places, most explicitly in Rcf. |(|. After this phase, the universe goes to the weak- field 
behavior, where our model behaves as usual. 

The main departure of our model from the usual one is in the region of large curvature scalar, where \B\ is large 



compared to a 3 . To find out the behavior of the solutions in this limit, we compare the leading terms in Eq. (7.17) 
and see that B is proportional to B 4 , where the coefficient of proportionality is positive and slowly varying. Thus, we 
find approximately B{a) w a 2 / 3 for small a. This implies the existence of a big-bang singularity, but with a different 
behavior: From Eq. ( [7.7[) we obtain 



ds 2 



da 2 

fl 2/3 



dx — dy 



dz< 



which corresponds in synchronized time to the behavior 



ds 2 = dt 2 - t 6/5 (dx 2 + dy 2 + dz 2 ), 



(7.18) 



(7.19) 



this being a good approximation to the exact metric for small t, differing from the usual big-bang behavior in almost 
all other models. Further details of our model will be presented elsewhere. 



C. The cosmological singularity 



Here we present the argument with the singularity behaviour mentioned at the end of section II: In our model, 
differently from Einstein's theory, the divergence of the curvature does not necessarily imply the divergence of any 



part of the energy- momentum-tensor. Let us concentrate on the trace. The r.h.s. of Eq. (2.11) reads 



R + 2l 4 R 3 

(i + ; 4 i? 2 ) 3/2 



8ttG 



- 3D 



(i + ; 4 i? 2 ) 



2\3/2 



±oo 



and this expression must be equal to the trace T of the energy-momentum tensor. In Einstein's theory, R 
necessarily implies T — > ±oo, whereas here, T may remain finite even if R — > oo. 

Detailed numerical calculations would support this qualitative picture, however, we postpone such calculations until 
we have a more strictly physically motivated form of the Lagrangian. 
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VIII. DISCUSSION 



We have argued that the gravitational action A has a decreasing dependence on R for increasing \R\. Such a 
behavior is expected from the spacetime stiffness caused by the vacuum fluctuations of all quantum fields in the 
universe. 

Our model does not have the tachyonic disease of R + R 2 models studies by Stelle and others fjli] . 

Since our model has an action which interpolates between Einstein's action and a pure cosmological term, it promises 
to have interesting observable consequences which may explain some of the experimental cosmological data. 

The heat-kernel expansion of the effective action in a curved background is closely related to the Seeley-Gilkey 
coefficients psj| , and for higher loop expansion also higher powers of curvature appear: To get the n-loop approximation 
one has to add terms until ~ R n+1 , a behavior which also happens in the string effective action ^o). So, if one cuts 
this procedure at a certain value of n, one gets always as leading term for high curvature values a term like ~ R n+1 . 
However, the n-loop approximations need not converge to the correct result if one simply takes n — > oo in the 
n-loop-result. In fact, what we have used in the present paper is such an example: 

c " °° k^gTTTFw °° TsSg + £ *' R " +1 M 

with some real constants fci,, where 



fc=i 



h = -^c < 8 - 2 > 

but the Taylor expansion on the right hand side diverges for R > l~ 2 . So, the Taylor expansion is useful for small 
i?-values only, and for large values R we need a correct analytical continuation. 

Prigogine et al. have proposed in Eq. (18) of Ref. |^TJ a model where the effective gravitational constant depends 
on the Hubble parameter of a Friedmann model. Though this ansatz depends on the special 3+1-decomposition 
of spacetime, it shares some similarities with the model discussed here. More recent developments how to find a 
well-founded gravitational action from considering quantum effects can been found in [^2f and p3|] . 

Quite recently, see for instance |p6| , accelerated expansion models of the universe have been discussed and compared 
with new observations. We postpone the comparison of our model with these observations to later work. 
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APPENDIX A: NEWTONIAN LIMIT IN A NONFLAT BACKGROUND 

The Newtonian limit of a theory of gravity is defined as follows: It is the weak-field slow-motion limit for fields 
whose energy-momentum tensor is dominated by its zero-zero component in comoving time. Usually, the limit is 
formed in a flat background, and sometimes, this is assumed to be a necessary assumption. This is, however, not true, 
and we show here briefly how to calculate the Newtonian limit in a nonflat background, Moreover, our approach is 
different from what is usually called Newtonian cosmology. To have a concrete example, we take the background as 
a de Sitter spacetime. 
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The slow-motion assumption allows us to work with static spacetime and the matter, assuming the energy- 
momentum tensor to be 



Tij = P 8^1 (Al) 

where p is the energy density, and time is assumed to be synchronized. The de Sitter spacetime in its static form can 
be given as 

ds 2 = -(1 - kr 2 )dt 2 + -^-r + rW, (A2) 
1 — kr z 

where x° = t, x 1 = r, x 2 = \, x 3 = 8 and dQ 2 — d\ 2 + sin 2 xdO 2 is the metric of the 2-sphere. In this Appendix, we 

have changed the signature of the metric from (H ), which is usual in cosmology, to ( — h ++), which leads to 

the standard definition of the Laplacian. 

The parameter k characterizes the following physical situations: For k = 0, we have the usual flat background. By 
setting k — Qwe can therefore compare the results with the well-known ones. The case k > corresponds to a positive 
cosmological constant A. In the calculations, we must observe that the time coordinate t fails to be a synchronized 
for k ^ 0, but it is obvious from the context how to obtain the synchronized time from it. 

In the coordinates (A2), there is a horizon at r = ro = So, our approach makes sense in the interval < r < tq. 
However, ro shall be quite large in comparison with the system under consideration, so that we do not meet a problem 
here. 

Now, the following ansatz seems appropriate: 

ds 2 = -(1 - fcr 2 )(l - 2^)dt 2 + (y^2 + rW ) C 1 + 2 ^)> ( A3 ) 

where if and ip depend on the spatial coordinates only. The weak-field assumption allows us to make linearization 
with respect to tp and An extended matter configuration can be obtained by superposition of point particles, so we 
only need to solve the problem for a <5-source at r = 0. This one is spherically symmetric, so we may assume ip — <p(r) 
and tp = tp(r) in Eq. (A3). For the metric components we get: 

.9oo = -(l-fcr 2 )(l-2^), gn = 1*1% 522 = r 2 (l + 2^), <7 33 = .922 ■ sin 2 X - (A4) 
The inverted components are up to linear order in tp and tp: 

9°° = ~~%~~THi> 3 11 = (l-fcr 2 )(l-2^), g 22 = 1 —^, 5 33 = 5 22 sin" 2 x , (A5) 
1 — kr r 

which gives the Christoffcl symbols 

*~» r -Tr5* < A6 > 

rj = (1 - kr 2 ) [-kr + 2kr{p + ip) - tp'(l - kr 2 )] , (A7) 

r 12 = r 3 3 = v' + -, (A9) 

r 

Tl 2 = -r(l - kr 2 ) -iP'r 2 {l - kr 2 ), (A10) 

r^sin 2 ^, (All) 

rl 2 = cot X , (A12) 

r 2 3 = -sinxcosx, (A13) 
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and the Ricci tensor reads 

Roo = -3k(l - kr 2 ) - tp"(l - kr 2 ) 2 - —(1 - kr 2 ) + 6k(tp + tp)(l - kr 2 ) + kr(l - kr 2 )(5p' - tp'), (A14) 

r 

R n = -2tp" + tp" - ~ r i>' + + Y ^- i W - V), (A15) 

R 22 = 3kr 2 - tp"r 2 (l - kr 2 ) - tp'(2r - Akr 3 ) + {tp' - tp')(r - kr 3 ), (A16) 

i? 33 = R 22 ■ sin 2 x- (A17) 

Before we discuss these equations, we consider two obvious limits: 

For k = 0, we see that i?oo = — <p" — 2tp' jr — —Ay, leading to the usual Newtonian limit Aip — —AnGp. 

For ip = ip — we get for the Ricci tensor: 

#g = R{ = Rl = Rl = 3k, (A18) 
and thus the de Sitter spacetime with R — 12k for k > 0. 

Returning to the general case we have 

— = 6k - I2k^: + [tp" - 2ip")(l - kr 2 ) + -tp' - 5krip' - -tp' + 7krtp' (A19) 
2 r r 



and then 



R 4 
R°- — = -3k + 6fc?/> + 2tp"(l - kr 2 ) - 6fcr^' + -ip'. (A20) 



The other components have a similar structure and can be calculated easily from the above equations. The first term 
of the r.h.s., — 3k, will be compensated by the A-term. The usual gauging to ip — > and tp — > as r — > oo is no 
more possible because for r > ro our approximation is no more valid. As an alternative gauge we add such constant 
values to tp and tp that they are approximately zero in the region under consideration. So we may disregard the term 
6kip. All remaining terms with k can be obtained from those without k by multiplying with factors of the type 1 + e 
where e ~ kr 2 , k = l/r 2 , with ro being of the order of magnitude of the world radius. In a first approximation, this 
gives only a small correction to the gravitational constant. In a second approximation, there are deviations from the 
1 /r-behavior. 

An analogous discussion for the Lagrangian R + l 2 R 2 tells us that in a range where I <C r <C ro, the potential 
behaves like (1 — c\e~ r ' l )/r, as in flat space. 



APPENDIX B: THE ABSENCE OF GHOSTS AND TACHYONS 



Here we show in more details what has been stated after Eq. (3.3). In the conformally transformed picture with a 
scalar field, the absence of tachyons (i.e., particles with wrong sign in front of the potential term) becomes clear from 
the form of the potential. For checking ghosts (i.e., particles with wrong sign in front of the kinetic term) we have to 
go a little more into the details: In Stelle |Q the particle content of fourth order gravity with terms up to quadratic 
order has been determined, and the existence/absence of ghosts and tachyons has been given in dependence on the 



free constants of the theory. In the first of Refs. [|4|] , the analogous calculation as in 27 has been done for a term 
R 3 added to the Einstein-Hilbert-Lagragian. Let us give here the argument for general n > 3: If R n is in C, then 
the term R n ~ l and its derivatives are in the corresponding expression after variational derivative with respect to the 
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metric. In the result, all terms represent products of at least n — 1 small quantities; because of n > 3 these are always 
at least two factors; thus, they all vanish in the linearization about the Minkowski space-time. 

Now, one might be tempted to require the analogous linearization properties for a Friedmann-Robertson- Walker 
background. However, linearization around other than fiat space-times is not at all a trivial task, see pq], even for 
Einstein's theory: For the closed Friedmann model, Einstein's theory is linearization unstable, for spatially flat models 
it is stable, and for the open Friedmann model the result is - contrary to other claims in the older literature - not yet 
known. We face the further problem that linearization around the de Sitter space-time is complicated to determine, 
because the same geometry can be locally represented as a spatially flat as well as a closed Friedmann model. So, we 
leave the question of linearization stability with non-flat background of our model unanswered. 

Another type of reasoning was given quite recently: In pjJ the possibility has been discussed that the contributions 
to the Lagrangian coming of gravitons on the one hand and of gravitinos on the other may cancel each other to avoid 
the ghost problem. 
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